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We study the “Higgs” amplitude mode in the relativistic quantum 0{N) model in two space dimen¬ 
sions. Using the nonperturbative renormalization group and the Blaizot-Mendez-Galain-Wschebor 
approximation (which we generalize to compute 4-point correlation functions), we compute the 0{N) 
invariant scalar susceptibility at zero temperature in the vicinity of the quantum critical point. In 
the ordered phase, we find a well-defined Higgs resonance for N — 2 and N = 3 and determine 
its universal properties. No resonance is found for > 4. In the disordered phase, the spectral 
function exhibits a threshold behavior with no Higgs-like peak. We also show that for N = 2 the 
Higgs mode manifests itself as a very broad peak in the longitudinal susceptibility in spite of the 
infrared divergence of the latter. We compare our findings with results from quantum Monte Carlo 
simulations and e = 4 — (d -|- 1) expansion near d = 3. 
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I. INTRODUCTION 

Relativistic quantum field theories with 0{N) sym¬ 
metry arise in the low-energy description of many 
condensed-matter systems: quantum antiferromagnets, 
superconductors, Bose-Einstein condensates in optical 
lattices, etc. In the ordered phase, where the 0{N) sym¬ 
metry is spontaneously broken, mean-field theory pre¬ 
dicts N — 1 gapless Goldstone modes corresponding to 
fluctuations of the direction of the A^-component quan¬ 
tum field, and a gapped amplitude “Higgs” mode (see, 
for instance, Ref. 1). 

The existence of the Higgs mode near the quantum 
critical point (QCP) separating the ordered and disor¬ 
dered phases, when fluctuations are taken into account 
beyond mean-field theory, has been a subject of debate. 
Does the Higgs mode exist as a resonance-like feature or 
is it overdamped due to its coupling to the Goldstone 
modes? Space dimensionality d plays a crucial role. In 
three dimensions, the QGP corresponds to a Gaussian 
fixed point of the renormalization group and interactions 
are suppressed at low energies; the Higgs resonance be¬ 
comes sharper and sharper as the QCP is approached. 
This has been beautifully confirmed in the quantum an- 
tiferromagnet TlCuCla."’^ By contrast, for 2 < d-l-I < 4, 
the QCP corresponds to the Wilson-Fisher fixed point 
and interactions are strong at low energies; the existence 
of a well-defined Higgs resonance is not guaranteed. Fur¬ 
thermore the visibility of the Higgs mode strongly de¬ 
pends on the symmetry of the probe.”* For 2 < d-|- I < 4, 
emission of Goldstone bosons leads to an infrared di¬ 
vergence in the longitudinal susceptibility,'’^* which is 
the standard correlation function to probe the amplitude 
mode, thus making the observation of the Higgs reso¬ 
nance very difficult. The 0(A^)-invariant scalar suscep¬ 
tibility (i.e. the correlation function of the square of the 
order parameter field) has a spectral weight which van¬ 
ishes at low energies and is a much better candidate.* 


The Higgs mode in a two-dimensional system has been 
observed in a Bose gas in an optical lattice in the vicin¬ 
ity of the superfluid-Mott-insulator transition’* and in a 
disordered superconductor close to the superconductor- 
insulator transition. *** 

The Higgs mode near a two-dimensional relativistic 
QCP has been studied with various techniques: large- 
N expansion,** quantum Monte Carlo simulations,*^”** 
nonperturbative renormalization group (NPRG),** and 
e = 4 — (d -I- 1) expansion about d = 3.*'’ For N = 2, 
these studies have conclusively shown the existence of 
a Higgs resonance in the ordered phase which persists 
arbitrary close to the QCP. However, besides quantita¬ 
tive issues, such as the precise value of the mass of the 
Higgs mode, some basic qualitative questions remain: is 
there a Higgs-like resonance also in the disordered phase? 
Does the Higgs resonance exist for larger values of N, e.g. 
TV = 3 and = 4? 

In this paper we address these issues using the NPRG 
approach.**”’^*’ More specifically, we use the Blaizot- 
Mendez-Galain-Wschebor (BMW) approach, an approx¬ 
imation scheme to the exact RG equations which allows 
one to obtain the momentum and frequency dependence 
of correlation functions.^*”^* The BMW approximation 
is also the starting point of the NPRG study reported 
in Ref. 16. In that work, however, additional (uncon¬ 
trolled) approximations were made in order to simplify 
the numerical solution of the NPRG-BMW equations. 

The main results of our analysis, which is restricted 
to zero temperature, are the following, i) In the or¬ 
dered phase, we hnd a well-defined Higgs resonance for 
Af = 2 in agreement with previous works. *^”**" ii) We 
show that a Higgs resonance is also present for A^ = 3 
but not for A^ = 4. This differs from previous NPRG 
analysis** but agrees with Monte Carlo (MG) results. *“*’*5 
iii) In the disordered phase we find that the spectral 
function of the scalar susceptibility does not exhibit a 
Higgs-like resonance peak above the absorption thresh- 
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old. This contradicts previous NPRG results^® and some 
QMC simulations.'^^ On the other hand, it is corrobo¬ 
rated by a separate MC analysis,'® which found no con¬ 
clusive evidence for such a resonance, and is supported 
by the results of the e = 4 — (d -|- 1) expansion, iv) Al¬ 
though the longitudinal susceptibility diverges as 1/w at 
low energies due to its coupling to the Goldstone modes, 
for N = 2 the Higgs mode manifests itself in the spectral 
function as a very broad peak. 

The outline of the paper is as follows. In Sec. II we 
present the NPRG approach to the quantum 0{N) model 
at zero temperature. In that limit, the two-dimensional 
quantum model is equivalent to the three-dimensional 
classical model. We show that the BMW approximation 
can be extended to compute 4-point correlation functions 
such as the scalar susceptibility. We also recall the ex¬ 
pected scaling behavior of the longitudinal, transverse 
and scalar susceptibilities near the QCP. In Sec. Ill we 
show that the BMW approach becomes exact in the limit 
N ^ oo. The numerical solution of the RG equations for 
N 1 and = 2, 3,4, • • • is discussed in Sec. IV. Both 
the longitudinal and scalar susceptibilities are considered. 
We compute the spectral functions and the correspond¬ 
ing universal scaling functions, as well as the mass of the 
Higgs mode and the stiffness. In Sec. V we provide a 
summary and conclusion. 

II. NPRG APPROACH 

The two-dimensional quantum 0{N) model is defined 
by the (Euclidean) action 

+ + (I) 

where we use the shorthand notation 

x=(r,r), J=J dr J (fr. (2) 

ip{x) is an V-component real field, r G [0,/3] an imagi¬ 
nary time and r a two-dimensional coordinate (/I = 1/T 
and we set h = = !)■ tq and Uq are temperature- 

independent coupling constants and c is the (bare) ve¬ 
locity of the (f field. The factor 1/N in Eq. (1) is in¬ 
troduced to obtain a meaningful limit V —>■ oo (with uq 
fixed). The model is regularized by an ultraviolet cutoff 
A. In order to maintain the Lorentz invariance of the ac¬ 
tion (I) at zero temperature, it is natural to implement 
a cutoff on both momenta and frequencies. 

The phase diagram of the quantum 0{N) model with 
V > 2 is well known. At zero temperature, there is 
a quantum phase transition between a disordered phase 
(^0 > ?' 0 c) and an ordered phase (ro < rpc) where the 
0{N) symmetry of the action (1) is spontaneously bro¬ 
ken (uo and c are considered as fixed parameters). The 


QGP at ro = rpc is in the universality class of the three- 
dimensional classical 0{N) model with a dynamical crit¬ 
ical exponent z = 1 (this value follows from Lorentz in¬ 
variance); the phase transition is governed by the three- 
dimensional Wilson-Eisher fixed point. At finite tem¬ 
peratures, the system is always disordered for V > 2, in 
agreement with the Mermin-Wagner theorem. For N = 2 
and ro < roc, there is a finite-temperature Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition^®^^® and the 
system exhibits algebraic order at low temperatures. The 
BKT transition temperature line Tbkt terminates at the 
QGP ro = roc- 

In the following, we consider only the zero temperature 
limit where the two-dimensional quantum 0{N) model 
is equivalent to the three-dimensional classical 0{N) 
model. For convenience, we set the velocity c equal to one 
so that the action (1) takes the usual form of the classical 
0{N) model with x a three-dimensional space variable. 
Having in mind the two-dimensional quantum 0{N) 
model, we shall refer to the critical point of the three- 
dimensional classical 0{N) as the QCP. In Fourier space, 
a correlation function xiPx^Py^Pz) computed in the clas¬ 
sical model should be identified with the correlation func¬ 
tion x{Px,Py, *w) in the quantum model (with w a bosonic 
Matsubara frequency) and yields the retarded dynami¬ 
cal correlation function x^iPx,Py,^) = xiPx^Py^i^ 
uj + i0+) after analytical continuation zw —>■ w -I- z0+. In 
the following, we shall often use the notation X'^(w) for 
X^{Px =Py = 0, w) and x"(w) = Im[x^(a;)]. 

A. Scale-dependent effective action r* 

To implement the NPRG approach, we add to the ac¬ 
tion an infrared regulator term, 

= i^(/?i(-q)i?fc(q)(/3i(q), (3) 

q.i 

such that fluctuations are smoothly taken into account as 
the momentum scale k varies from the ultraviolet cutoff 
A down to The cutoff function in (3) is defined by 

Rfe(q) = Zkq^r , r{y) = (4) 

where a is a constant of order one (we choose a = 
2.25) and Zk a field renormalization factor defined in 
Sec. HD.-'® Thus the ASk term suppresses fluctuations 
with momenta |q| ^ A: but leaves unaffected those with 
111 ^ A;. We consider the ^-dependent partition function 

Zk[3,h]= (5) 

in the presence of external sources J and h. The mean 
value of the field is defined by 

J,/i] = = (V3(x)). (6) 
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The bilinear source h will allow us to compute the scalar 
susceptibility (but will eventually be set to zero). 

The central quantity in the NPRG approach is the 
scale-dependent effective action 


Tk[(t),h] = - In Zk[3,h]+ / 3-(l>-ASk[(p], (7) 

J X 

defined as a modified Legendre transform of — In Zk [3,h], 
wrt the linear source J, which includes the subtraction 
of In Eq. (7), J(x) = h] is obtained by 

inverting Eq. (6). Assuming that for k = A the fluctu¬ 
ations are completely suppressed by the ASk term, we 
have 


rA[cf>,h] = S[ct,]- [ hcf>‘^. ( 8 ) 

J X 

On the other hand, the effective action of the original 
luodel (1) is given by Tfc^o since Rk=o vanishes. The 
NPRG approach aims at determining Pfe^o from Pa using 
Wetterich’s equation*^' 

d^Tu[<t>M = ^TrjRfc (Pf’°)[(/.,h] + , (9) 


where Rk = d^Rk and denotes the second- 

order functional derivative of Pfe [(p, h] with respect to cp. 
In Fourier space, the trace involves a sum over momenta 
as well as the 0{N) index of the p field. 

Most of the physically relevant information encoded in 
Pfc can be obtained either from the effective potential or 
the one-particle-irreducible (IPI) vertices. The effective 
potential is nothing but Pfc evaluated in a constant, i.e. 
uniform, field configuration (up to a volume factor E), 


Ukip) 


y^kiP] 


cf)—const 


( 10 ) 


Because of the 0{N) symmetry of Pfc, Uk is a function 
of the 0{N) invariant p = 0^/2. We denote by po.fc 
the value of p at the minimum of the effective potential. 
Spontaneous symmetry breaking of the 0{N) symmetry 
is characterized by a nonvanishing expectation value of 
the field y, i.e. limfc_>.o po,k = Po > 0. In this case, the 
effective potential Uk=o{p) is constant (and minimum) 
for 0 < p < Po in agreement with the convexity of the 
Legendre transform Pfc=o [</’]• At the QGP (ro = roc), 
Po = 0 but po,fc > 0 for all nonzero values of k. 

The IPI vertices are defined by 


_ s-+^rk[p,h] _ 

^Pii (xi) • • • 5pi^ (x„)i5/i(yi) • • • Sh{ym)' 

The correlation functions evaluated for h = 0 and in a 
uniform field configuration are determined by the vertices 






(n,m) 

k,{ij} 


[{xi},{yj };</>, h] 


h=0 


( 12 ) 


In particular, the propagator Gk,ij{p,p) = 
{pi{p)p’j{—p))h=o in a uniform field is obtained 
from the matrix equation 

Gk{p,P) = {r^k’°\p,P) + Rky\ (13) 

where p[.^’°^(p, <^) = (p, —p, ^). The 0{N) sym¬ 

metry allows us to write 

(P:^) = SijPA,kiP, P) + PiPj^B.kiP, p), (14) 

which yields the longitudinal (L) and transverse (T) parts 
of the propagator, 

Gfc.L(p,p) = [PA.fc(p,p) + 2pPB.fe(p,p) +Rfc(p)]"\ 
Gfe.T(p,p) = [PA.fc(p,p)+ 77fe(p)]"^ (15) 

Due to rotation invariance in space, 2-point vertices and 
correlation functions in (14) and (15) are functions of 
P=|P|- 

Important information can be obtained from the lon¬ 
gitudinal and transverse susceptibilities that we define 
by 

Xa(p) = Gfc=o.a(p, Po,k=o) (a = L, T). (16) 

In the disordered phase (po = 0), xl(p) = Xt(p) = 
Xl,t(p)- The single-particle excitation gap A is obtained 
from the vanishing of the spectral function %((rp(p,a;) = 
Im[xL rp(p, w)] for |w| < A. The excitation gap manifests 
itself as a sharp peak in t(‘^) = Xl t(p = (®®® 

Sec. IV B 2). In the ordered phase, the stiffness Ps is 
defined by*^ 

Xt(p) = for p ^ 0. (17) 

PsP^ 

For two systems located symmetrically wrt the QGP (i.e. 
corresponding to the same value of |ro — rod), fo ^be 
ordered phase (with stiffness ps) and the other in the dis¬ 
ordered phase (with excitation gap A), the ratio Ps/A is 
a universal number which depends only on N. This al¬ 
lows us to use A as the characteristic energy scale in both 
the disordered and ordered phases (in the latter case, A 
is defined as the excitation gap at the point located sym¬ 
metrically wrt the QGP).^' Note that A and Ps vanish 
as |ro — rod'^ as we approach the QGP. 

In the universal regime near the QGP (scaling limit), ^ 

Xa(p) = (I-) , 

(18) 

X"(w) = Im[xf (w)] = , 

where rj is the anomalous dimension of the if field at 
the QGP. $ct,± and are universal scaling func¬ 

tions and Za^± a nonuniversal constant with dimension 
of (length)''. The index +/— refers to the disordered 
and ordered phases, respectively. At the QGP (A = 0), 
Xa(p) -- P""' and xd(‘^) -- Since XlH and 

Xx(<^) are odd in uj we shall only consider the case w > 0 
in the following. 
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B. Scalar susceptibility 


We now consider the scalar susceptibility 

Xsiy - y') = (¥’(y)V(y')^) - (‘^(y)^)(‘^(y')^) 

S'^\nZ[3,h] 


6h{y)dh{y') 


(19) 


J—h—0 


where Z[3, h] = Zk=o[3, h]. Using (7), we can express 
as a functional derivative of the effective action F = F k=o. 


ShiyysHy') 


( 20 ) 


where the order parameter (p[h] is defined by 


sr[cf>,h] 

(50 (x) 


= 0 . 

4>=<j>[h] 


( 21 ) 


In Eq. (20), 5/5h{y) is a total derivative which acts both 
on (f)[h] and the explicit h-dependence of the functional 
T[cf>,h]. Using (21), one finds 


S^r[(f)[h],h] 

Sh{y)6h{y') 


F(O’2)[y,y';0[/,],/j] 


+ 


^( 1 . 1 ) 


x,y;0[/i],/i] 


S4>i [x; h] 
Sh{y') 


( 22 ) 


(we use Einstein’s convention for summation over re¬ 
peated indices). To compute Scj)[h]/Sh, we take the func¬ 
tional derivative of Eq. (21) wrt h, which gives 


<50i[x;h] 

Sh{y) 


- J ^[x,x';0[h],/i] 

X Tf''^\x',y;^[h],h], 


(23) 


where ^ is the propagator. From (20,22,23), we 

finally deduce 


X(5(y-y')= - r^°’^^(y,y',0)-f f Ff’^^(x,y,0) 

J x,x' 

X (x, x', 0)F^.^’^^ (x', y', 0) (24) 

(with 0 = 0[h = 0]) or, in Fourier space, 

Xs(p) = -r(°’^)(p, 0) 

+ Ff’1) (p, 0)Fg’°)-i (p, 0)Ff’1) (p, 0). (25) 

The last term corresponds to the part of the scalar sus¬ 
ceptibility which is not IPI. Equation (25) is shown dia- 
grammatically in Fig. 1. 

Since 0 transforms like a vector under 0{N) rotations 
while h is a scalar. 


(a) Xs ~ -f -w*. 


{b) = 





FIG. 1. Diagrammatic representation of the scalar suscepti¬ 
bility (25) (a). A vertex is represented by a black dot 

with n solid lines and m wavy lines (b). Solid lines connecting 
vertices stand for the propagator G = 


where / and 7 are functions of p = |p| and the 0{N) 
invariant p. To determine the scalar susceptibility in 
the NPRG approach, we must therefore consider the k- 
dependent vertices and or, equivalently, the 

fc-dependent functions /fc(p,p) and 7 fc(p,p). 

In the universal regime near the QCP,^^ 

x«(p) = s±+yi±A3-2/-|.,_±(|-), 

(27) 

X"(a;) = Im[xf (a.)] = (|) , 

where and are universal scaling functions and 
A±,B± nonuniversal constants. At the QCP (A = 0), 
Xs(p) - X(((0) ~ and x"(w) - gj^j^e 

Xs (w) is an odd function of uj, we shall only consider the 
case w > 0 in the following. 


C. BMW approximation 

In this section, we first review the BMW approxi¬ 
mation for the solution of the RG equation (9) when 
h = We then show how this approximation can be 

extended to the calculation of the scalar susceptibility. 


1. Effective potential Uk and 2-point vertex 

Equation (9) cannot be solved exactly. However, it 
can be used to derive flow equations for the effective po¬ 
tential and the IPI vertices. The effective potential is 
determined by 

dkUkip) = \ J ^fe(q)[G'fe,L(q, p) + {N - l)Gfc,T(q, p)] 

(28) 

with initial condition Ua{p) = rop-\- {uo/6N)p'^. We use 
the notation 


rf’^^(p, 0 ) = 0 */(p,p), 

F(°>2)(p, 0) = 7(p,p), 


( 26 ) 


Solving (28) requires to know the propagator Gk = 
+ Rk)~^ and therefore the 2-point vertex in 
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a uniform field. The latter satisfies the equation (Fig. 2) 


(P’ H ^kGk,i,^^ (q, (t>) 


(p> -p> q> -q> (P’ “P “ 

^ (30) 


ds.o) 


X Gk,^3^^ (p + q, (l>)^k’Zn (“P’ p + q> -q> ‘^) 


where 

dkGk^iii2 (Q; 4^) — (Q)G^/c,zii3 (Qi 4^)Gk,i^i2 (Q: 0) ■ 

(31) 

More generally, the IPI vertices satisfy an infinite hier¬ 
archy of equations since involves and 

The BMW approximation^allows us to 
close this inhnite hierarchy of equations. It is based on 
the observation that the i?fe(q) term in (31) restricts the 
integration over the loop momentum in (30) to small val¬ 
ues |q| < k whereas the regulator term ASk ensures that 
the vertices r[,"’°^(pi • • -p^,^) are regular functions of 
the momenta when |pi| ^ k. This enables us to set 
q = 0 in the vertices and in (30). Noting 

that in a uniform field, 



FIG. 2. Diagrammatic representation of the RG equa¬ 
tions (30,39,40). Signs and symmetry factors are not shown. 
The diagrammatic representation of the vertices is the 

same as in Fig. 1 and the cross stands for dkRk- 


equations is obtained by writing the self-energies as^^ 


^aAp^p) = P^^A,fe(P,P), 

AB,fc(P,p) = P^TB,fc(p,p), 


(35) 


p(3,0) 

k,111213 


(P: 


-p, 0 , 0 ) 


p(4,0) 

^ k ,11121314 


(P, 


-p, 0 , 0 , 0 ) 


1 -p, 0 ) 

Av 90,3 

1 92 r£°)jp,-p, 0 ) 

^ 90,3 90 ,, 

(32) 


we obtain a closed equation for the 2-point vertex. 


and solve for Ya^r and Ys^k with initial conditions 

YaMp^p) = ^b,a(p,p) = 0. 


2. Scalar susceptibility 


, ,.,92rg|°)(p,0) 


9fcri^^^(P:^) = X / [^kGk,il^2{<^,4>)] 


d4>iid(t)i, 


9r£g(p» 

90i2 

d^tfhpA) 


- / [dtGk,1112 A, (p)] 


X Gfc,,3„ (p-1-q, 0 )— ‘ ' 1 (33) 


The vertices and F^’^^ in a constant field satisfy 
the RG equations 


9 fer^°’^^ (p, (p)=Y^ dkGk ,1112 (q, 0 ) 


P’ “P’ “ ^kSi3 “P “ P’ 


which must be solved together with (28). 

Equation (33) yields two coupled equations for 
rA,fc(p,p) and TB,k{p,p) [see Eq. (14)]. Eollowing 
Refs. 21-23, we introduce self-energies AA,k and AB,k 
defined by 


X Gk,i 3 u (p + q, (P + “P’ ’ 


(36) 


and 


rA,fc(p, p) = P^ + A^,fc(p, p) + uAp), 

rB,fc(p,p) = AB,fe(p,p) -h u'Ap), 

with AA,fe(p = 0,p) = AbAp = 0>p) = 0- The effec¬ 
tive potential Uk is obtained from the (exact) RG equa¬ 
tion (28) while the self-energies are deduced from the 
(approximate) RG equation (33). The final form of the 


9feriy ^ (p, 0 ) = X! 9 fcGfc,, 3,3 (q, 0 ) 

q 


X 


IpfS,!) 

2 k,ii2ii 


(p, q, -q, -p, 0 ) - PkfX (P’ “P “ 
(p + q,0)ri^)^q(p + q,-q,-p,0) • 


q, 0 ) 
(37) 
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The BMW approximation amounts to setting q = 0 in 
the vertices in (36) and (37). Using 


riYUp,o,-p,^) = 




rSiL3(P’'^’0,-p,^) = 


1 5r[^^)i^(p,-p» 


'Jv 

1 52r^°’^^(p,-p,</)) 


V d<j3i^d(j)i^ 
1 


, (38) 


one then obtains closed equations for and 

4r. (p.« = 2 / 




5ri^J(p,(/)) 


X G/c,i3iJp + q, (t>) 


and 


^fcriY'(P:^) = ^ / 5fcGfc.iii2(q,<A) 


d4>Z2 

d^PkfiP^(f>) 


(39) 


d(t)i^d4>i 


/"a^ ^^^^£S(P’'^) 

Jq ^4^12 

^ / , ^^^^^!M^(p,(/)) 

X Gfc,i 3 i Jp + q, (/))- ’— -. (40) 

Equations (39) and (40) are shown in Fig. 2. They lead to 
RG equations for the functions fk and defined in (26) 
with initial conditions /a(p,p) = —2 and 7a(p,p) = 0. 


D. Explicit form of RG equations 


the k, p and p dependence of the functions (and use the 
notation p = |p|). The 2-point vertex is determined by 
the flow equations 

^ [Js'^ip^Yi + W'f + Jj^ip^YB + W')^ 

- + 2pY'l) 

-\ir[{N-l)Y'^ + 2YB\ (43) 

and 

= ^{(^ - i)[jr{p'YB + w'f - irw^] 

+ J3^^[p^(iA + 2 Ub -f 2pyi^) -t 3W' -f 2pW"f 

- + 2pW''f - jY^{p^y^ + W'f 

+ J3LTp^'2 _ + w'f + if^W^] 

- \ir{N - l)Y{, - + 2py") 

^ ^{ir - I2 ^)Yb. (44) 

2p 

Note that the rhs of the last equation is well defined for 
p —>■ 0 and p -A 0.^^ Equations (42-44) were previously 
derived in Refs. 21-23. 

/ and 7 satisfy the equations 

dtf = Jt^if + 2pf')[pfYi + 2 Yb + 2pY{,) + 31U' 

+ 2pW"] + (TV - fjriiP^YB + W') - I^^pf" 

-lm2^ + {N-l)ir] (45) 

and 


da = (TV - i)jrf + + M'f 

- ^ [l2^{f + 2P7") + {N- l)irf]- (46) 


In this section, we give the explicit forms of the flow 
equations in the BMW approximation. We define the 
threshold functions 

Jk!liPW)= [ [9tT?fc(q)]G”^(q,p)Gfe.p(p-tq,p), 

J q 

Ik^niP) = Jk!'niP = 0w), (41) 

where a,l3 = L,T and t = ln(fc/A) is a (negative) RG 
“time”. 

The derivative Wk{p) = U'j.{p) of the effective potential 
satisfies the equation 

dtW = i [II^' + (TV - 1 )/TT'], ( 42 ) 

where the primes denote a p derivative. Here and in 
the following, to alleviate the notations, we do not write 


The QGP manifests itself as a fixed point of the RG 
equations provided we use dimensionless flow equations 
where all quantities are expressed in units of the running 
scale k. We therefore use the following dimensionless 
variables (with d = 3) 

P = k~^p, p = K~^Zkk'^~‘^p, (47) 

and functions ®^ 

Uk{p)=K-^k-'^Uk{p), 

Wk{p) = Z-^k-^Wk{p), 
l + 4(4,fc(p,p) = Z-fl + YA,k{p,p)], 

YB,k{P^ P) = KdZ^^k'^~‘^YB,kiP,p), 
fk{p,p) = fk{p,p), 
lk{p,p) = K-^Zlk‘^-‘^-fk{p,p), 
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where = 2‘^“^(i7r‘^/^r((i/2) is a constant originating 
from angular integrals which is introduced here for con¬ 
venience. The field renormalization factor is defined 
by 


Zh. = 




dp^ 


= 1 -h yA,fe(p = o,po.fc), 


P=0,p=po,fe 


(49) 

where po.fe corresponds to the minimum of the effective 
potential Uk{p)- Equation (49) implies that 


renormalized in the large-iV limit. The other equations 
read 

dtWu = ^ir'. 
dtTB,k=Njrri^k - 

N 

dtfk=NjrhTB,k - -jiriL 
dak = 


r[,^;^^(p,po.fc) - ^feP^ + bEfc(po.fe) for p^O (50) 

so that in the ordered phase the stiffness is simply defined 
by Ps,k = ‘^ZkPo,k [Eq. (17)]. If we use (50) to estimate 
the gap A in the disordered phase (where po,k=o = 0), 
we find 


( Wk=o{0) 

V Zk=o 


1/2 


(51) 


The numerical solution of the flow equations show that 
this expression is in very good agreement (with an er¬ 
ror smaller than l%o) with the exact determination of 
the gap obtained from the peak in the spectral function 
xIt{oj) (Sec. IVB2). 

We also define a (running) anomalous dimension 


rjk =-kdklnZk- (52) 

At criticality (rg = rgc), the anomalous dimension is 
obtained from rj = limfe_>.o pk. The normalization con¬ 
dition (49) can be expressed as YA,k{p = ^,Po,k) = 0 
which provides us with an equation for pu- This equa¬ 
tion takes a simpler form if we choose p = 0 in the defi¬ 
nition (49) of This choice is however not appropri¬ 
ate in the ordered phase where we are eventually inter¬ 
ested in quantities defined at the nonzero order parame¬ 
ter Pox, which corresponds to a diverging dimensionless 
field Pox ~ Po,k/k. 


III. LARGE-A LIMIT 

In this section, following Ref. 21, we show that the 
BMW equations become exact and can be solved analyt¬ 
ically when N ^ oo. In this limit, T^)^] is of order N 
and the field 4> is of order y/W (see Appendix A). This 
implies that pk is 0{N), Wk, Tax fk are 0(1) (as 
well as the threshold functions and whereas 

Tbx is 0{1/N). It follows that 


To obtain the equation for Tbx, have used = 
and = —21T^/j’'^ when Aax = 0- To 

solve these equations, we set W = Wkip) and use the 
variables (fc, W) instead of (fc, p).-' This is done introduc¬ 
ing the function gkiW) = p and using (/(.(IT) = l/IT^(p), 


dt9k{W) = -^^d,Wk{p) = ^ir 

_ N r 9 t 7 ?fc(q) 

“ 2 4 h^ + W + RkiqW 


(55) 


Since k and W are independent variables, this equation 
can be rewritten as 


dtgk{W) 


f 1 

2 ‘4q^ + IT + Rfe(q)’ 


(56) 


where both sides are total derivatives, and we obtain 


gk{W)-gA{W) 


^[( _^_ 

2 J^W + W + Rk{q) 


q^ + W + Rxiq) 


(57) 


Eor i?A(q) —t 00 , using gxiW) = {iN/uo){W — To) we re¬ 
produce the known result in the large-A limit [Eq. (A8)]. 
Eor i?A(q) < oo, we obtain an apparent difference with 
the exact result, which is explained by the fact that Ea)*/*] 
is not given by £'[(/>]. This does not matter when one is 
interested in universal properties in the vicinity of the 
QCP: the microscopic physics can be directly parame¬ 
terized by Ta)*/*] (which we can choose to coincide with 
S')*/)]). In that case, however, it is important to keep the 
term with Ra in (57) for comparison with the numerical 
solution of the flow equations in the large-A limit. 

To calculate Tbx, 

dtTBx\w= dtTBx\p+^Bx9tgk = NJ^Tlx ( 58 ) 

(the prime denotes a p derivative), i.e. 

dtTB\\^=-Njr. (59) 


dtYAx 


A 


tTTy' 

-'2 ^AX' 


(53) 


The rhs is again a total derivative. Using T b,a = mo/ 3A, 
we obtain 


Since Ya,a = 0, we deduce Yax = Aax ~ 9: the mo¬ 
mentum dependence of the transverse propagator is not 


TBxiP,P) 


1 

^ +'f [nfe(p,p) -fiA(p,p)] 


(60) 











where 
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nfc(p,/o) 


1 


/q q2 + Vhfe(p) + i?fc(q) 

1 


(P + q)^ + Wk{p) + i?fc(p + q) ■ 


(61) 


We recover the large-A^ result derived in Appendix A 
when i?A(q) —t 00 and IIa —>■ 0. 

A similar procedure is used to calculate fk and 7^. 
From 


dtfk\^ = NjrrB,kfk, 

dak\^ = Njrfl 

and the initial conditions /a = — 2, 7 a = 0, we obtain 


fk{p,p) 

7fc(P,p) 


6N 

Uo 


^B,k{p,p) 



2 

^B,kiP,p), 


(63) 


in agreement (when i?A(q) —t 00) with the re¬ 
sults (A15,A18) obtained in the standard large-A^ ap¬ 
proach. 

The numerical solution of the flow equations in the 
large-A^ limit is discussed in Sec. IV A. 


IV. HIGGS AND LONGITUDINAL SPECTRAL 
FUNCTIONS 

We solve numerically the flow equations with uq = 200 
and A = 1. We use a (p,p) grid of 80 x 200 points with 

^ — P Pniax; 0 ^ P ^ Pmax5 Pmax — 8 and Pmax — 

6N. The flow equations are integrated using an explicit 
Euler method with At = —10“"^ {t = ln(A:/A)). We use 
Simpson’s rule to compute momentum integrals (with an 
upper cutoff ^max = 4) and finite-difference evaluation for 
derivatives. We have verified the stability of our results 
wrt to the various parameters used (number of points in 
the {p, p) grid, upper cutoff in momentum integrals, etc.). 

The ordered phase requires some care. First, since 
po,fc converges towards a nonzero value po for fc —)■ 0, 
Po,fc ~ Po/^ diverges and one cannot work with a fixed p 
grid. To circumvent this difficulty, we start the flow with 
a fixed p grid but switch to a fixed p grid (0 < p < Pmax) 
once the flow leaves the critical regime to reach the or¬ 
dered regime, which occurs for k of the order of the in¬ 
verse of the Josephson length ^ (roc — ro)~''■ Second, 
since Wkipo,k) = 0, Wk{p) is negative for p < po.fe and 
Gfc(p = 0,p = 0) = (i?fc(p = 0)-|-VFfc(0))“^ becomes very 
large. This behavior is associated to the approach to con¬ 
vexity of the effective potential Uk{p)^^ While Gfe(p,p) 
should remain strictly positive, we find that it develops 
a pole for small k in the BMW approximation, which 
leads to a divergence of the threshold functions (41) 
and a numerical instability.^^ To cure this problem, we 



■NPRG 
exact i?A < oo 
exact i?A = oo 


0 ^ 
10“^ 


10 " 


10 " 


p/A 


FIG. 3. (Color online) Scalar snsceptibility Xa)?) at the QCP 
for N = 1000 (solid line) compared to the exact large-A solu¬ 
tion with Ra < oo (circles) and Ra = oo (squares). Here and 
in the following figures we use arbitrary units for the scalar 
and longitudinal susceptibilities. 


eliminate from the p grid the points for which 114 (p) 
is smaller than —2.^® The grid [pmin.fe; Pmax] becomes k 
dependent.When Pmin.fc becomes nearly equal to po,fe, 
the flow cannot be continued anymore since the mini¬ 
mum of the effective potential must remain in the range 
[Pmin.fc, Pmax] for physical quantities (defined for p = po.fc) 
to be determined. However, this procedure allows us to 
reach small values of k (typically fcmin — 0.05A) for which 
physical quantities have nearly converged to their k = 0 
values (the only exception is the very low-energy behav¬ 
ior p, jw] <C A of correlation functions, see below).'"*® A 
more precise estimate of the fc = 0 values can be obtained 
using an extrapolation of the type a + bk'^. 

A momentum dependent function F(j>) = Fk=o{p), 
such as a two-point vertex or the scalar susceptibility, 
is obtained from the approximation F{p) ~ Fk^p/p^^^{p) 
where k = p/pmax is the smallest value of k for which 
the dimensionless momentum p = p/k is still in the grid 
[OjPmaxj- We have verified, by increasing Pmaxj that the 
flow of F’fc(p) for k < p/pmax is negligible. This is due to 
the fact that, in the cases we are considering here, p acts 
as an effective infrared cutoff, while only momenta of the 
order of k or smaller contribute to the flow, so that the 
flow of the function (p) effectively stops when k p. 
To obtain the retarded dynamical function F^{uj), we 
compute F{p) for M momentum values pi (Z = 1 • • • M) 
with typically M in the range 50 — 100. One then con¬ 
structs a M-point Fade approximant Fp{p) which coin¬ 
cides with F{p) for all pi’s, and F^{uj) is approximated 
by Fp{—ioj + e) (we take e/A < 10“'*).®’’ Note that in 
the ordered phase, we cannot determine F{p) for values 
of p below fcminPmax = S^min — 0.4A. This prevents us 
to determine F^{uj) for w < 0.4A. On the other hand, 
in the disordered phase, the spectral functions we are in¬ 
terested in vanish for w < A or a; < 2A and there is no 
need to compute F{p) and F^{uj) for p ot uj smaller than 
0.4A. 
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FIG. 4. (Color online) Spectral functions ’Xsiu) and Xl(‘^) 
in the ordered and disordered phases for N = 1000 (solid 
line), compared to the exact large-solution (symbols). In 
the disordered phase, the exact solution for Xl(cu) = Xt{^) ~ 
(5(a; — A) is not shown. 


TABLE I. Critical exponent i/ obtained in the NPRG ap¬ 
proach, from either Wk{Q) (Sec. IVBl) or Xs (Sec. IVB3), 
compared to Monte Carlo (MC) simulations. 


A 

from 144(0) 

from Xs 

MC 

2 

0.673 

0.674 

0.6717(1) [38] 

3 

0.714 

0.722 

0.7112(5) [39] 

4 

0.754 

0.766 

0.749(2) [40] 

5 

0.787 

0.804 


6 

0.816 

0.835 


8 

0.860 

0.879 


10 

0.893 

0.906 


100 

0.990 

0.992 


1000 

0.999 

0.999 



TABLE 11. Same as Table I but for the anomalous dimension 
t; (defined by ry = limfe_>o rjk when ro = roc)- 


A 

NPRG 

MC 

2 

0.0423 

0.0381(2) [38] 

3 

0.0411 

0.0375(5) [39] 

4 

0.0386 

0.0365(10) [40] 

5 

0.0354 


6 

0.0321 


8 

0.0264 


10 

0.0220 


100 

0.00233 


1000 

0.000233 



A. Large-A limit 


As a check of our procedure we first discuss the nu¬ 
merical solution of the flow equations in the large-A limit 
where comparison with exact results is possible (Sec. III). 
Figure 3 shows the scalar susceptibility Xs(p) obtained 
for N = 1000 at the QCP. Except for momenta near the 
cutoff A, we obtain a very good agreement with the ex¬ 
act solution (63) in the limit A —> oo taking into account 
the finite value of i?A(p)- For sufficiently small p, when 
nfe=o(p) becomes very large, the scalar susceptibility be¬ 
comes independent of the initial value Rjy of the cutoff 
function. In any case, for universal properties, the value 
of i?A does not matter. 

The spectral functions Xs(^) ^-nd Xl('^) 3-^® shown in 
Fig. 4 for A = 1000, in both the ordered and disordered 
phases, in the universal regime near the QCP. Again, the 
agreement with the exact results (including nonuniversal 
prefactors) in the limit A —>■ oo is very good. This val¬ 
idates our procedure to compute the momentum depen¬ 
dence of correlation functions as well as the Pade method 
to obtain the spectral functions. 


B. A < 00 

In the following, we discuss the NPRG results obtained 
for finite A, in particular A = 2 and A = 3. 


1. QCP 

We first solve the equations to determine rgc and the 
critical exponents v and p. The anomalous dimension rj 
is directly obtained from limfe_>o?7fc when rg = rgc. The 
correlation-length exponent ly is deduced from the be¬ 
havior of Wk{0) — IT*(0) -I- at very long time \t\ 

(since the condition rg = rgc is never exactly fulfilled the 
RG trajectories will always eventually flow away from the 
fixed point with an escape rate given by 1 /v). Our results 
agree both with previous NPRG-BMW calculations^^ 
and Monte Carlo estimates (see Tables I and II). 

Figure 5 shows Xl,t(p) and Xs(p) at criticality for 
A = 2 and A = 3. In the universal regime p po, where 
Pg is the inverse of the Ginzburg length ^g ~ 247r/Mg, we 
find Xl,t(p) \/p^~'^ where the value of the anomalous 
dimension, rj ~ 0.0423 for A = 2 and p ~ 0.0411 for 
A = 3, agrees with the estimate obtained from the run¬ 
ning anomalous dimension (Table II) . As for the scalar 
susceptibility, we find Xs(p) ~ with 9 ~ 0.0345 for 
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FIG. 5. (Color online) |Xs(p) — Xs(0)|/p®~^^‘" and 
Xl,t(p)/p~^^'’ at the QCP for N = 2 and N = 3. The nor¬ 
malization is chosen to have a ratio equal to one for p —>■ 0. 


N = 2 and 9 ~ 0.230 for N = 3. If we use the expected 
relation 9 — 3 — 211/ [Eqs. (27)], we obtain i/ ~ 0.674 for 
N = 2 and v ~ 0.722 for TV = 3, in very good agreement 
with our previous estimates of v based on the behavior 
of Mffc(O) in the close vicinity of the fixed point (Table I). 


2. Disordered phase 


Figures 6 and 7 show xl,t(p) and Xs(p) and their 
spectral functions in the disordered phase for N = 2 and 
N = 3. The various curves, obtained for different values 
of ro — roc, show a data collapse in agreement with the 
scaling forms (18,27) expected in the critical regime. The 
excitation gap A, deduced from the peak in the spectral 
function Xlt('^)j ia very good agreement with the 
approximate expression (51). 

The spectral function x"(w) of the scalar susceptibility 
vanishes for jw] < 2A. Contrary to previous conclusions 
based on QMC and we find that 'x's{uj) 

rises smoothly above the threshold at w = 2A with no 
sign of a local maximum for lo > 2A. The authors of 
Ref. 15 argued that in spite of the maximum observed 
above the threshold in their MC simulations, there is in¬ 
clusive evidence for a resonance at finite frequency in the 
disordered phase (the peak carries a small spectral weight 
and its position is not very robust). We also note that 
no resonance is obtained in the 4 — (d -|- 1) expansion. 




p/A 




FIG. 6. (Color online) Xl,t(p) and Xl,t(<^) in the disordered 
phase for N = 2 and N = 3. The solid line and the symbols 
correspond to different values of ro — roc. 




p/A p/A 



FIG. 7. (Color online) Same as Fig. 6 but for the scalar 
susceptibility Xs. 


3. Ordered phase 


In Fig. 8, we show the derivative Wk{p) of the effec¬ 
tive potential for various values of k and N = 2. As ex¬ 
plained above, for small k we must use a fc-dependent grid 
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p 


FIG. 8. (Color online) Derivative Wk{p) of the effective po¬ 
tential, shown in the range [pmin.fc, Pmax], for various values of 
k (from bottom to top: ln(fc/A) = —11.5, —12.25, —13 and 
— 15.25) and N = 2. The inset shows the corresponding flow 
of po.fc vs ln(fc/A). 


TABLE 111. Universal ratio ps/{NA) obtained from the 
NPRG in the BMW approximation (NPRG BMW). Also 
shown are the previous NPRG results obtained from a deriva¬ 
tive expansion of the effective action"*^ (NPRG DE) and from 
Monte Carlo simulations’ ’ (MC). The exact result in the limit 
N —>■ 00 is I/Itt ~ 0.0796. 



ui/A ui/A 


FIG. 9. (Color online) Xs(p) ^^nd Xs (‘^) in Ih® ordered phase 
for N = 2 and N — 3. The solid line and the symbols corre¬ 
spond to different values of ro — roc- 


N _1000 10 8 6 4 3 2 

NPRG BMW 0.0796 0.0803 0.0829 0.0903 0.111 0.137 0.193 
NPRG DE^*’ 0.0838 0.085 0.086 0.089 0.096 0.106 0.132 
MC’"’ 0.114 0.220 


[Pmin.fc, Pmax] to get rid of the smallest p values for which 
the propagator is not positive. For k = fcmin — 0.05A, 
Pmin,k IS very close to po,fc and we cannot continue the 
flow. In Fig. 8 we also show the behavior of po.fe and its 
convergence towards its fc = 0 value. The extrapolated 
value at A: = 0 differs from the value at fcniin by less than 
1 %. 

Table 111 shows the universal ratio Ps/{NA) where p^ 
and A are computed for the same distance |ro—rod to tbe 
QCP (see the discussion at the end of Sec. 11 A). For small 
values of N we find significant deviations wrt previous 
NPRG results."” The value for A = 2 is now much closer 
to the MC estimate of Ref. 15 and the agreement is also 
satisfactory for A = 3. For A = 1000, we recover the 
large-A result. 

In Fig. 9 we show Xs (p) and x"(w) in the ordered phase 
for A = 2 and A = 3. Again, we observe data collapse 
in agreement with the scaling forms (27). For A = 2, 
we find a well-defined Higgs resonance whose position 
u! = mu and full width at half-maximum vanishes as the 
QCP is approached. For mn u> pc^ recover 
the critical scaling x!sipj) ^ Up to a multiplica¬ 

tive factor which depends on the nonuniversal factor A- 
[Eq. (27)] the shape of the resonance, given by the univer- 




I 

<1 




u/A 


cj/A 


FIG. 10. (Golor online) Log-log scale plot of Xsd^) In the 
ordered phase for A = 2 and A = 3, showing the asymptotic 
behavior Xs (u)) ~ ca® at low energies. 


TABLE IV. Universal ratio mujA obtained from the NPRG 
in the BMW approximation. Also shown are previous NPRG 
results’® as well as results obtained from (Q)MG”’’'’ and e 
expansion.’ ' 


A 

3 

2 

NPRG BMW 

2.7 

2.2 

NPRG’® 


2.4 

MG’-* 

2.2(3) 2.1(3) 

QMG’’ 


3.3(8) 

e expansion’’ 

1.64 

1.67 
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FIG. 11. (Color online) Spectral function Xs{^) for various 
values of N in the ordered phase. 


sal scaling function is in very good agreement with 

the MC result of Refs. 14 and 15. The Higgs resonance 
is still visible, although less pronounced, for iV = 3. This 
observation disagrees with previous NPRG results^® but 
agrees with MC simulations of Ref. 15. The universal ra¬ 
tio itih/A, shown in Table IV, is compatible with MC es¬ 
timates of Refs. 14 and 15. Since in the ordered phase we 
must stop the flow at a finite value k^in, we cannot cal¬ 
culate reliably the spectral function x"(cu) for frequencies 
u! < fcmiij. Although for fcniin ^ w < mn, our results are 
compatible with x'si^) (see Fig. 10), the low-energy 
regime w <C A where the spectral function is completely 
determined by the Coldstone modes is difficult to access. 
In Fig. 11 we show Xsi^^) Af = 2, 3,4, 5,10,100. Only 
for V = 2 and N = 3 does a Higgs resonance exist. 

Finally we show the longitudinal susceptibility xl(p) 
and its spectral function Xl(<^) Fig. 12 for iV = 2 and 
N = 3. For p —>• 0, the longitudinal susceptibility Xl(p) 
diverges as 1/p as expected for a two-dimensional sys¬ 
tem (Fig. 13).'*^ This effect is a consequence of the cou¬ 
pling of the longitudinal mode to the Coldstone modes''^”® 
and prevents the observation of a well-defined Higgs reso¬ 
nance in Xl(^)-^ Nevertheless a broad peak, presumably 
due to the Higgs mode, can be seen for w ~ mn when 
N = 2 (Fig. 12).'*^ For N = 3, the peak has disappeared 
but a faint structure can still be seen. 


V. CONCLUSION 


We have studied the scalar and longitudinal suscepti¬ 
bilities in the quantum 0{N) model using the NPRC. 
Comparison with QMC simulationsand e = 4 — (d-l- 
1) expansion^^ allows us to identify robust properties of 
the Higgs mode: i) In the ordered phase, there is a well- 
defined Higgs resonance for N = 2 and N = 3. The 
spectral function x/i^) ti^s been determined both from 
QMC and NPRC but the precise value of the the mass of 
the Higgs mode is not precisely known (Table IV). If we 
take the difference between NPRC and MC simulations 




p/A 




j/A 


FIG. 12. (Color online) Xl(p) and Xl(^) in the ordered phase 
for N — 2 and N — 3. The solid line and the symbols corre¬ 
spond to different values of ro — roc- 


10 ° 

r<i 

<1 

10-1 


lO"! 10° 
oj/A 

FIG. 13. (Color online) Log-log scale plot of Xl(^) in the 
ordered phase for A = 2 and A = 3, showing the asymptotic 
behavior Xl('^) ~ 1/^ at low energies. 



of Refs. 14 and 15 as an estimate of the error, then the 
ratio niH /A is known within 5% for A = 2 and 20% for 
N = 3. ii) In the disordered phase, there is no Higgs-like 
peak in Xsi^) above the absorption threshold. There are 
two other important properties obtained from the NPRC 
that have not been studied with MC or other methods 
so far: iii) The Higgs resonance is suppressed for A > 4. 
iv) For A = 2 the Higgs mode manifests itself in the lon¬ 
gitudinal spectral function Xl('^) by a very broad peak. 

From a more technical point of view, we have shown 
that the BMW approximation^!-^® allows one to compute 
the momentum dependence of correlation functions, in¬ 
cluding 4-point functions such as the scalar susceptibility. 
We have also shown that the difficulties arising from the 
non-positivity of the propagator in the ordered phase can 
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be overcome by using a fc-dependent grid which does not 
include small values of the order parameter p. 
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Appendix A: Vertices in the large-A limit 

In this section, we determine the effective potential 
Wk and the vertices and in the large- 

N limit using the standard approach where the partition 
function Zk [h] = [h, J = 0] is obtained from a saddle- 

point calculation.We first introduce the field p = 
and a Lagrange multiplier A, 


Zu[h]= j V[p,p,\]eyipiy- j pf + p + ^ - p) c^(x)i?fc(x - x') • ¥’(x') 




^(2/i + iX- rof - i [(Vv’)^ H- iXp^] 


f ¥3(x)Rfc(x-x') • ¥>(x')j. (Al) 


Then we split the field p into a a field and an {N — l)-component field tt. Integrating over the tt field, we obtain the 
action 


Sk[cr,X,h] = f 

J X 


-‘^{2h + iX-roY + \ [(^cr)^ -p lAcr^] 


+ \J f7-(x)Rfe(x-x')cr(x')-b^^^^y^Trlngj,i[A], (A2) 


where 


Let us first consider the effective potential for h = 0, 


5 ^^[x,x';A] = [-V^-|-*A(x)]J(x-x')-|-Rfe(x-x') (A3) 

is the inverse propagator of the field in the fluctuating 
A field. In the limit N —>■ oo, the action becomes pro¬ 
portional to N (if one rescales the a field, a —)■ '/Na)\ 
the saddle point approximation becomes exact for the 
partition function Zk[h] and the Legendre transform of 
the free energy coincides with the action Sk ■ This im¬ 
plies that the scale-dependent effective action (7), defined 
as the Legendre transform including the subtraction of 
AS'fe[cr], is simply equal to Sk{<J, A, h] — AS'fc[(T]: 


rk[a,X,h] = y 1“^ (2^+ 

+ 2 I + (A4) 


(we use V — 1 ~ A for large N). We can eliminate the 
Lagrange multiplier field using 


Srk[cr,X,h] 

SX(x) 


= 0 , 


A=Afc[CT,/i] 


(A5) 


to obtain a scale-dependent effective action P^ [cr, h] = 
L k [o’, Xk[cr, h ], h] which is a functional of a and h. Vk{cr,h] 
is the starting point to compute the vertices in the 

large-A limit. 


Uk{p) = - X —{iXk - ro)^ -I- iXkP 

2uo 

A r 

+ Y J + -Rfc(q)], (A6) 

where we use the notation p = cj^ 12. The value of 
Xk = Xk{p) is obtained from dUk/dXk = 0, which fol¬ 
lows from (A5), i.e. 

- {iXk-ro)+p+— / ., , „ ; Y = 0- (A7) 

uq 2 jq -l- iXk + i?fc(q) 

We deduce that 


Wk{p) = iXk 


= ro 


uop 

3A 


uq r 1 

6 J^q^ + Wk{p) + Rkiqy 


(AS) 


which is the known result in the limit A —)■ 

The vertex can be obtained from Tk[cr, X] setting 

( 2 ) 

h = 0. In Fourier space, the 2-point vertex T), can be 

( 2 ) ( 2 ) 

written as a 2 x 2 matrix with components F^ 

^ k,X(T and Inverting this matrix, one obtains the 

^‘2^\ _ 

longitudinal propagator as (Ty )cro-(p) = [rA./c(p) + 
2prB,fc(p)]~^ from which we deduce (see Eq. (55) in 
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Ref. 8) 


rA,fc(p,/o) = 


rB.fc(p,p) = 


3 ^ ^TT / 

— + ^nfc(p,p) 

Uo 2 


-1 


(A9) 


The momentum dependence of the transverse 2-point 
vertex rfc,T(p,p) = P^ + Wk{p) remains the bare one. 

Let us now consider the vertex (x, x', p), 




5'^Tk[(J, Xk[cr, h],h] 


Sh{x)Sh{x.') 


, (AlO) 


<j{z)=a,h=0 


where cr = and d/Sh{x) denotes a total derivative 
(Sec. lie 2). Using (A5), we obtain 




S‘^Tk[a, Xk[(J,h],h] 


(5h(x)(5/i(x') 
d'^Tk[a,X,h] (5Afc[y; cr, h] 


a(z)=a,h=0 


Sh{x.)6X{y) Sh{x.') 


(T{z)—a,X—Xk [<7,0] ,/i=0 

(All) 


and in turn 


ri°’^\x,x',p) = - ^^S{x-x') 

Uq 

6N i5Afc[x; cr, h] 


12N 


— i 


Uq Sh(x') 


(A12) 


CT(z) = tT,/t = 0 


To determine the last term of (A12), we take the func¬ 
tional derivative 6/6h{x') of Eq. (A5), which gives 


0= -i—S{x-x') 
Uo 


[ rSA(x,y,p) 

-'V 


(5Afe[y; cr, h] 


where 


,p) = 


Sh{x') 

S^rk[a,X,h] 


, (A13) 


<j{z)=a,h=0 


sX{x)6X{x') 


(T{z)—a^X—Xk [cr,0],/i.=0 

(A14) 

From (A12,A13), we finally deduce 


r 

where 


12iV , n(2) 


Uo 


k (p,p) = —:— + ( — ) r^;,(p,p) , (A15) 


Uo J 


rSA(p>p) = ^ + ynfe(p,p), 

Uo ^ 


nfc(p,p)= / 9k{ci,p)gk{p + q,p) 

J a 


(A16) 


and 5fc(q,p) = [q^ -fi Wk{p) + i?fe(q)] ^ (we use iXk = 
Wkip)). ^ 

Following a similar approach for 


riY^(x,x',p) = 


one finds 


5‘^Tk[<7, Afc[cr, h],h] 


6a{x)Sh{x') 


^k,i\p^p) = -^V^^k,i\iP^P) (A 18 ) 

Uo 


, (A17) 


<j{z)=a,h=0 
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